#ifndef SMOOTHING_H
#define SMOOTHING_H

#include "symphonie.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_bspline.h>
#include <gsl/gsl_multifit.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_statistics.h>

  #define N        6


/* number of fit coefficients */
#define NCOEFFS  8

/* nbreak = ncoeffs + 2 - k = ncoeffs - 2 since k = 4 */
#define NBREAK   (NCOEFFS - 2)

int smooth_Bspline(QList<float> X,  QList<float> Y, float val);

#endif // SMOOTHING_H
